Hückel theory is a simple and powerful method for predicting the molecular orbital and the energy of conjugated molecules. However, the presence of nitrogen atoms in aza aromatic molecules alters the Coulomb and resonance integrals owing to the difference in electronegativity between nitrogen and carbon atoms. In this study, we focus on acridine and phenazine. Further correction is implemented based on the ring current model, thus revealing the change in resonance integral for the carbon-carbon bond along the bridge of the molecule. The Hamiltonian of the π-electron system in the Hückel method is solved using the HuLiS software. Various geometry-based aromaticity indices are used to obtain the aromaticity indices of the two non-equivalent rings. For further evaluation, the results for bond lengths are used to calculate the associated bond energy. Considering the carbon-hydrogen (CH) bonds, the total molecular energy is compared with the experimental heats of formation for a number of benzenoid hydrocarbons and aza aromatics, in addition to the two studied molecules. Finally, the correlation between the nitrogen atom on the aromaticity index and the ring energy content is evaluated to determine to which extent the Hückel model agrees with previous experimental and advanced computational studies. However study phenazine carboncarbon molecule πelectron π electron software geometrybased geometry nonequivalent non equivalent rings evaluation carbonhydrogen hydrogen CH (CH bonds aromatics Finally studies

Excessive degeneracy in the ground state of the π-electron energy levels of anthracene was removed using the Hückel method by correcting the Coulomb integral of the four central carbon atoms. A further correction to the resonance integral was proposed based on the ring-current model, which describes the π-ring current flow along the molecule’s perimeter, by introducing a new parameter for the bridge carbon atoms expressed as a fraction of the resonance energy. The solution to the Hamiltonian was obtained based on the symmetry group theory, which provides the advantage of solving the determinant matrix or secular equations in a particular irreducible representation with a relatively small matrix dimension. The results were further analyzed to determine the bond order and bond length. The values of the harmonic oscillator model of aromaticity, GEO, and EN indices of aromaticity and their ratio for the central and outer rings of benzene were used to evaluate the validity of the correction. Applying the same method to phenanthrene as a topological analog of anthracene allows for an interesting comparison of the stabilities of the two molecules.